x1 = normrnd(5,1,50,1);
x2 = normrnd(6,1,50,1);
x = [x1 x2];
boxplot(x,1,'g+',1,0)
